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Abstract. We discuss a thin film evolution equation for a wetting evaporating liquid on a 
smooth solid substrate. The model is valid for slowly evaporating small sessile droplets when 
thermal effects are insignificant, while wettability and capillarity play a major role. The model 
is first employed to study steady evaporating drops that are fed locally through the substrate. 
An asymptotic analysis focuses on the precursor film and the transition region towards the bulk 
drop and a numerical continuation of steady drops determines their fully non-linear profiles. 

Following this, we study the time evolution of freely evaporating drops without influx for 
several initial drop shapes. As a result we find that drops initially spread if their initial contact 
angle is larger than the apparent contact angle of large steady evaporating drops with influx. 
Otherwise they recede right from the beginning. 



1. Introduction 

Evaporation of thin liquid films and sessile droplets has attracted much atten- 
tion both as the way to probe the dynamics of the contact line Il30l [19 1 and 
as a route to create deposition patterns through sedimentation of solutes and 
suspensions HJ El El El. A number of studies concentrate on problems per- 
tinent to any evaporation process, which include mass and heat transfer and 
thermocapillarity 0QX For slowly evaporating small sessile droplets stud- 
ied in contemporary well-controlled experiments on smooth surfaces IT30l l5ll. 
thermal effects are, however, insignificant, while contact line dynamics and 
liquid-substrate interactions play a major role. It has been suggested that the 
relation between spreading and evaporation/condensation goes both ways, so 
that the latter may alleviate the notorious contact line singularity Il42ll26 1. For 
background on spreading see, e.g., Il4ll36l. 

A remarkable phenomenon observed in evaporating completely wetting 
liquids is the formation of a dynamic meniscus with a finite contact angle 
J8l|5l|32l. The standard approach to computing the form of a spreading and 
evaporating drop and the resulting dynamic contact angle (3l Q31 123 is based 
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on the lubrication approximation with the singularity at the contact line al- 
leviated by slip. The evaporation rate is treated in three alternative ways. 
One possibility, realized in the presence of a temperature difference between 
the substrate and vapor phase Q, is the evaporation rate determined by the 
balance of latent heat and heat flux. The evaporation rate is then uniform 
in the limit of small Biot numbers, but diverges near the contact line in the 
opposite limit. Both limits can also be realized under isothermal conditions. 
The evaporation rate is uniform (as long as the layer thickness remains outside 
the range of intermolecular forces) when evaporation is controlled either by 
phase transition kinetics at the interface or by diffusion through a boundary 
layer of constant thickness in a stirred vessel. If evaporation is diffusionally 
controlled with no stirring, the evaporation rate increases towards the contact 
line; an analytical solution based on interfacial equilibrium with no flux onto 
the unwetted substrate yields the flux diverging on the contact line Il9l l27ll . 
Another approach accounts for the influence of the thermal conductivity of 
the substrate and the dependence of the vapor saturation concentration on 
temperature lfT"3ll3Tll . 

The aim of the present paper is to present a simple isothermal thin film 
evolution equation with evaporation limited by phase transition kinetics (or 
boundary layer transfer, but not diffusion), that correctly describes the influ- 
ence of effective molecular interactions on evaporation in the case of com- 
plete wetting. This is achieved by taking into account the dependence of the 
saturated vapor pressure on the disjoining pressure and curvature in the way 
it has been done in studies of the dynamics of evaporating films ll24l l20l l25l 
but not in the cited studies of droplet spreading. This allows us to describe 
in a consistent way the transition from the bulk droplet to a precursor layer 
and eliminate singularities at the contact line. Note however, that our model 
may be obtained as the isothermal limit of the models in (II EH; i-e., letting 
the difference of substrate and ambient temperature, and the latent heat go to 
zero. It also corresponds to the limit of infinite thermal conductivity of the 
liquid. 

Following EUl, we consider a two-dimensional "fed" system that al- 
lows us to study steady states of evaporating droplets. These steady states are 
compared to droplet shapes resulting from a time evolution of an evaporating 
droplet (without influx). The comparison will be employed to discuss a pos- 
sible special role the steady state profiles play in the time evolution. A related 
approach is taken in Ref. where steady fronts of evaporating liquid on an 
incline are considered. 

The paper is structured as follows. The following section[2]introduces our 
model and discusses the scaling, whereas section[3]discusses the asymptotics 
in the precursor film. Section [4] discusses the properties of steady drops with 
influx as a function of the influx strength and of the single remaining dimen- 
sionless parameter. The time-evolution in the case without influx is analysed 
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in section [5] where we also compare the steady drop profiles in the case with 
influx to the time-dependent profiles in the non-fed case. Section [6] gives our 
conclusions. 



2. Basic Model and Scaling 

For simplicity, we restrict our attention to a two-dimensional system as sketched 
in Fig.[T] Conceptually, there exists no difference to the full three-dimensional 
system, we only expect the transport rates to change. 




Figure 1. Sketch of the two-dimensional geometry employed for investigating an evaporating 
droplet with localised influx q(x). 

Using the lubrication approximation, the evaporation dynamics for an isother- 
mal droplet of liquid on a porous substrate is captured by an evolution equa- 
tion for the film thickness profile h B51 12U1 1251 15511 

d t h = -d x jconvO) -ievap(^) + <?0), (1) 
iconv(x) = d x p, jevapix) = ft ~ /Xfj^ , 

P = ~ld xx h - 11(h). (2) 

The first and second term on the r.h.s. of Eq. ([T} are the divergence of the 
convective flux j conv (x) and the evaporative flux j eva , p (x) , which correspond, 
respectively, to the conserved and non-conserved part of the dynamics. The 
function q(x) is the influx through the (locally) porous substrate. The evap- 
orative flux is proportional to the difference between the chemical potential 
of the ambient vapour and the chemical potential in the liquid fi = p/p; 
p is pressure, (3 is an effective evaporation rate constant; and 7, p and 77 
are the surface tension, mass density and dynamic viscosity of the liquid, 
respectively. The pressure p contains the curvature pressure —jd xx h and 
the disjoining pressure n(/i) modelling wettability J6l [161 |23l \TE\ U|; the 
hydrostatic pressure is neglected as we focus on nano- and micro-droplets. 
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To model a droplet of completely wetting liquid, we employ a long-range 
stabilising van der Waals disjoining pressure II = — A/Qirh? with the Hamaker 
constant ^4 < (29l SSI [34l [20l [37] [381 . Note that other sign conventions are 
also common (cf. e.g. |6l|35l|4]]). The model ([TJ is related to various models 
in the literature: it may be obtained from the one in |[20l by adding an influx 
and replacing the disjoining pressure for a partially wetting liquid by one for 
a wetting liquid. The models in Q] [28l incorporate various thermal aspects 
that are here neglected by assuming that the latent heat is very small or/and 
the thermal conductivity is very large. The same applies to the steady state 
description in ETTl . Note that our model also corresponds to the one in ||28ll 
in the limit of zero superheat. 

A dimensionless form of Eqs. ([T]> and Q can be obtained by choosing 
the characteristic energy density of molecular interactions between the fluid 
and substrate k = |^4|/67r/ijj as the pressure scale and the equilibrium film 
thickness ho = \ A/Q-irp^ 1 ' 3 corresponding to the ambient vapour potential 
Ho as the scale of film thickness h (note that po < when a thick flat film 
evaporates). The horizontal coordinate x and time t can be scaled in several 
ways [25 ]. A short horizontal length scale 

jh l2 / 67r 7 f\ A \Y /6 \/7 ^ 




6-7T J |/9/io| 2//3 

is fixed by the balance between disjoining pressure and surface tension at 
the thickness of the wetting layer, and determines the extent of a region 
adjacent to the contact line where the interface may be strongly curved due 
to interaction with the substrate. The lubrication approximation remains for- 
mally applicable as long as / far exceeds ho- Note, however, that lubrication 
approximation often still predicts the qualitative behavior for many systems 
with larger contact angles ll23l [T8l . When considering the results obtained 
with models like Eqs. ([T]), one has always to keep in mind that even very large 
contact angles obtained in lubrication approximation (measured as slopes at 
the inflection point of the drop profiles) correspond to rather small angles in 
physical scaling. 

Another horizontal scale, applicable in the precursor layer, is determined 
by the balance of flow driven by the disjoining pressure gradient and evapo- 
ration: 




L=\HF = \ ~ 7T, (4) 



This scale is large when evaporation is slow. It is appropriate to choose L 
as the horizontal scale, assuming it to be of the same order of magnitude 
as the third available scale - the droplet size. The respective time scale is 
T = (L/ho) 2 rj/K, and the flux j CO nv is scaled by IiqL/T. 
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Retaining the same notation for the rescaled variables, we rewrite Eqs. ([I]) 
and (|2]) as 

d t h = -d x j com/ (x) + (ed xx h + ^ - 1^ + q(x), (5) 

h 3 ( 1 \ 

d x ed xx h + — r , (6) 



where the parameter 



3 x V hy ' 

(6vr) 2 / 3 7/3r/ 



1^12/3^4/31^11/3 



(V) 



denotes the scale ratio {1/ L) 2 . Note that it is proportional to the evaporation 
rate constant /3, but contains as well a weak dependence on the chemical 
potential /xo in the denominator. 

In the following we will study steady state droplets that are obtained for 
an influx q(x) localised at the centre of the drop (section [4]). Below, the 
steady profiles are compared to time simulations without influx for different 
initial profiles (section [5]>. First, however, we discuss the asymptotics in the 
precursor film. 



3. Asymptotics in the precursor film 

In the outer precursor region, the film is almost flat and surface tension can be 
neglected, i.e., we set e = in Eqs. ^ and ^ and assume q(x) = 0. In the 
linear regime the film thickness decays exponentially to its equilibrium value 
h = l: 

h — 1 ~ exp(— V3x). (8) 

Note the difference from the non-physical asymptotics h ~ x 1 / 4 in E71 
where the dependence of the evaporation equilibrium on the disjoining pres- 
sure was neglected. The latter profile corresponds to the well-known result of 
de Gennes [6] who failed to recognise it as an unstable solution. 

To discuss the nonlinear behaviour we reduce Eqs. ([5]) and (|6]) to the 
stationary equation 

dx 2 /i 3 

This equation is solved by using h as an independent variable, and y(h) = 
(ciln h/dx) 2 as a dependent variable. The transformed equation is 

m _ Ui - . do* 
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It is integrated with the boundary condition y(l) = to yield 



y(h) 



In h 



1 f 1 

si 1 "* 



The precursor film profile is obtained in an implicit form 



V2a 



H - - (1- e- 3H ) 
3 v ; 



-1/2 



dH, 



(11) 



(12) 



where H = Inh. This solution formally implies a very fast growth h ~ 
exp[(x — xq) 2 ] towards the bulk of the droplet. It becomes, however, inappli- 
cable as h grows, necessitating a modified scaling. One can see that the two 
terms in the r.h.s. of Eq. ([6]) become, up to logarithmic corrections, compara- 
ble at h ~ e -1 / 4 , which, though appreciably exceeding the thickness of the 
equilibrium wetting layer h = 1, may be still far below the height of the bulk 
droplet. As follows from Eq. ( [TT] ), the incline at this thickness level is, up to 
logarithmic corrections, h x = h^y(h) ~ e~ 4 / 4 in agreement with results 
by Morris ||2T1 l22l . This sheds light on the origin of a finite contact angle 
in an evaporating droplet. As we will see below in section [4] the numerically 
obtained dependence agrees well with the asymptotic result. 

The flux J from the droplet bulk into the precursor at a "transitional" 
location X corresponding to the thickness level h = e -1 / 4 £ is determined by 
the total evaporation rate from the precursor, which can be obtained directly 
from Eq. Q: 



J 



x 



h{x) c 



dx 



dlnh 

dx 



x=X 




a/4 



(13) 

The dependence both on e and on a precise choice of the level £ is very weak. 
The rest of evaporation goes at an almost constant rate from the bulk of a 
large droplet. 



4. Steady state droplets with influx 

For zero influx through the porous substrate (q(x) = 0) the only steady state 
solution is h = h^. However, for q(x) ^ steady droplets may exist with 
a volume determined by the dynamic equilibrium between the overall influx 
through the substrate and the overall evaporation flux. 

Here we use continuation techniques lTT2l [TUl fTTTl to numerically analyse 
the steady state solutions of Eqs. ([5]) and ([6]), i.e., we set dt h = and solve 
the resulting ordinary differential equation as a boundary value problem on 
a domain of size D with the boundary conditions (for a symmetrical drop) 
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d x h = d xxx h = at x = (drop centre). At x = D we employ either h = 1 
and c^/i = 0, or d x h = and d xx h = 0. If the domain is sufficiently large the 
results depend neither on the particular choice of D nor on the used version 
of boundary conditions at x = D. For details on the usage of continuation 
methods for thin film equations see, e.g., Refs. fiH . (T7\ and j39l where they 
have been employed to study sliding drops, chemically driven running drops 
and drops pinned by wettability defects, respectively. 
For the influx q(x) we use a normalised Gaussian 

,.2" 



<?o- 



OV7T 



cxp 



x- 



(7 



(14) 



with qo = J °° q{x)dx being the total influx through the substrate. If the 
droplet size is large as compared to the width a, the results do not depend on 
the particular choice of a. 
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Figure 2. (colour online) For the case of small droplets we give (a) droplet profiles and (b) 
evaporation flux in dependence of position for e = 1.0 and various total influxes go as given 
in the legend. Domain size is D — 10, and a — 0.1. 




0.1 0.2 0.3 0.4 0.5 0.2 0.4 0.6 0.8 1 



x/D x/D 
Figure 3. (colour online) For the case of small droplets we show (a) droplet profiles and (b) 
evaporation flux in dependence of position for go = 2.5 and various e as given in the legend. 
Domain size is D — 10, and a = 0.1. 

Figs. [2] and [3] show profiles of rather small (nano-)droplets (left panels) 
and the corresponding dependencies of the evaporative flux on position (right 
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panels). The results are given for various moderate values of the influx qo 
(Fig. [2]) and the length scale ratio e (Fig. [3). In all shown cases these droplets 
are not much higher than the wetting layer. For such small drops the behaviour 
is dominated by the influence of the disjoining pressure. In consequence, the 
evaporation decreases monotonically from the center of the drops towards 
the contact region. Interestingly, in all cases the droplet assumes a shape that 
does not allow for any condensation of liquid even in the contact line region 
where the Laplace pressure is negative. Note that there exists a one-to-one 
correspondence between the strength of the influx qo and droplet volume for 
fixed e. This implies that one may characterise the relative size of droplets 
either by volume or by influx qo. 

For extremely small drops (see, e.g., profile for qo = 0.5 in Fig. [2]) the 
disjoining pressure influence is stronger than the capillary pressure even at 
the drop centre. As a result, the absolute value of the evaporation flux jevapis 
smaller than one even at the centre of the drop. For slightly larger drops 
(see, e.g., profile for qo = 2.5 in Fig. [2]) the capillary pressure dominates 
the disjoining pressure at the drop centre and j ev ap is larger than one. With a 
further increase in drop size the influence of the capillary pressure diminishes 
and jovap eventually approaches unity everywhere with the exception of the 
contact line region (cf. Fig. [4]). 

Decreasing e mainly influences the height of the droplets while the width 
remains roughly constant [Fig. [3ja)]. This implies that the curvature at the 
drop centre and the apparent contact angle 9 app (defined as the maximal slope 
of the drop profile), increase with decreasing e. However, although curvature 
increases we find that the influence of capillarity on the evaporation flux 
decreases [Fig.j3|b)]. For e = 0.01, one has at the drop centre j eva ,p slightly 
above one. Furthermore, at the same e, j CV ap has already a small plateau at 
the drop centre, i.e., the flux is nearly constant at the value determined solely 
by the chemical potential. 

The influence of the source width a is marginal as long as it is sufficiently 
smaller than the droplet width. For moderately large width it has still no 
influence on the contact line region but has some influence on the center of 
the drop. Increasing, for instance, a from 0.1 to 1.0 at constant jo = 2.5 and 
e = l the drop volume goes up by about 5%. Decreasing a down to 0.001 has 
no visible influence on the droplet shape. 

The droplets discussed up to this point represent nano-droplets of heights 
normally below 500 nm. However, for much smaller e or much larger qo one 
is able to study micro-droplets with heights in the 10-100 /im range. Fig. ffl 
shows profiles of such drops and the local evaporative flux for e = 10~ . 
For such large drops the local evaporation is essentially constant for the 'bulk 
drop' and decreases monotonically in a confined contact region [Fig. |4jb)]. 
Panel (c) of Fig.|4]shows the logarithm of h — 1. By shifting the drops in the 
x-direction one can appreciate that the approach to h = 1 is universal and 
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Figure 4. (colour online) For the case of large drops we show for e = 10 -6 (a) drop profiles 
and (b) evaporation flux in dependence of position. Results are given for various total influxes 
jo and droplet volumes V (see legend). The thin dotted line in (a) gives for the largest drop 
the corresponding parabolic drop profile of identical maximal height and curvature at centre 
(corresponding to a spherical cap in lubrication approximation). Panel (c) shows log ft to 
indicate the universal behaviour near the contact line (drops shifted in x). The dotted line 
indicates the linear result h — 1 ~ exp(— y/3x) [Eq. jsjl]. Domain size is D = 50, and 
a = 0.1. 




Figure 5. Shown are (a) drop volume and (b) the apparent contact angle S app (defined as the 
maximal slope of the drop profile), in dependence of total influx qo for various length scale 
ratios e as given in the legend. The straight dotted [dashed] lines indicate linear [quadratic] 
dependencies, respectively. 
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well described by the linear relation h 
Eq. ([8]> in section [3]). 



1 ~ exp(— y/3x) derived above (see 




Figure 6. (colour online) Shown are (a) drop volume and (b) apparent contact angle in depen- 
dence of the length scale ratio e. The total influx go is constant for each line, respectively. The 
lines are characterised by the drop volume at e = 1.0 (see legends). 



When increasing the influx for fixed e the steady drops become larger in 
width and height [Fig. |2ja)]. This is indicated as well by the dependence of 
volume on influx [Fig. [5ja)]. The corresponding apparent contact angle is 
shown in Fig. [5jb). One clearly distinguishes a small-drop and large-drop 
regime with a crossover at about V = 1. In the small-drop regime volume 
and contact angle are both proportional to the influx. In the large-drop regime 
the contact angle approaches a constant (or increases with growing influx fol- 
lowing a power law with an exponent smaller than 1/5), whereas the volume 
depends quadratically on influx. The latter is easily explained noticing that 
the evaporative outflux for large drops is proportional to the surface "area" 
of the drop (negligible influence of Laplace and disjoining pressure). For a 
constant contact angle the area under the parabola depends quadratically on 
its arc length. For the influx to balance the outflux, the surface area has to 
grow proportionally with the influx, i.e., the volume increases quadratically 
with the influx. 

Inspecting Fig. [5]f urther, one notices that the overall behaviour is different 
for larger (e > 10 ) and smaller (e < 10~ 3 ) drops. In the former case the 
transition between the small-drop and large-drop regime is monotonic, i.e., 
the slopes of the curves in Fig.[5]change monotonically. In contrast, for small 
e < 10~ 3 in the transition range one may define a third region where the 
slopes of the V(qo) and # app (go) curves pass through a maximum. 

The tendency towards a constant contact angle for increasing volume can 
also be observed in Figs. [6ja) and (b) where we plot the drop volume and 
the apparent contact angle, respectively, as a function of the length scale 
ratio e for various fixed influxes for rather large drops. We find that for large 
drops, the volume as well as the contact angle decrease for increasing length 
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scale ratio e roughly as e -1 / 4 . This agrees with the asymptotic expression 
determined above in Section |3] 

For smaller drops, deviations from the power law are found at larger e. 
Interestingly, the dependence of the contact angle on e seems to approach a 
limiting curve for large drops. In the following, we employ the curve for the 
largest drops in Fig. [3] as an approximation to the asymptotic dependence of 
#app on e for infinitely large drops. 

5. Time evolution without influx 

Next, we study the time evolution of evaporating droplets without influx 
through the substrate, i.e., we simulate Eq. (|5]> with q{x) = 0. The domain 
size D and boundary conditions at x = and x = D correspond to the 
ones used in the steady state calculations in the previous section. We use 
three different initial profiles hi{x) = h(x,t = 0) of equal maximal height 
h m and volume V: (i) a parabola hi(x) = {h m — 1)(1 — x/x c ) 2 + 1 with 
x c = 3V/ {h m — 1) for < x < x c and hi{x) = 1 for x > x c \ (ii) a gaussian 
hi{x) = {h m — 1) exp((x/cr) 2 ) + 1 with a = IV j ^pn{h m — 1); and (iii) the 
steady-state solution of identical V and h m obtained in section|4j 




120 



Figure 7. Space- time plot of an evaporating droplet for e = 10~ 4 . The initial profile is a 
parabola on a precursor film of thickness h p = 1. It has a volume of V = 1000 and maximal 
height of // m ax = 140.6. The corresponding contact angle is 8[ n i — 26.3. The initial height 
corresponds to the one at V = 1000 for the steady state drops with influx for the corresponding 
e (obtained in section|4j. 

Fig [7] shows a space-time plot for a typical time evolution observed when 
using an initial parabolic drop profile that has the same height and volume as a 
fed drop obtained in sectional The case shown is for e = 10~ 4 . At early times 
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Volume Volume 
Figure 8. Trajectories in the phase plane spanned by the maximal drop height and drop volume 
for (a) e = 1(T 6 and (b) e = 1. Shown are curves resulting from (i) time evolutions for three 
different initial profiles of equal maximal height and volume (parabola, gaussian and steady 
state with influx), and (ii) calculations of steady state solutions with influx as obtained by 
continuation. 



the contact line region relaxes under the influence of the disjoining pressure, 
thereby decreasing the apparent contact angle. Subsequently, the width and 
height of the drop decrease monotonically until at about T = 100 the drop 
has vanished and only the stable precursor film remains. When starting (as 
in the present case) with the drop measures (volume and height) as obtained 
for the drop with influx, the evolution always looks similar. In particular, 
we have not found that the drop macroscopic ally spreads at the beginning 
(by "macroscopic" we mean a spreading that goes beyond the small local 
relaxation at the contact line). 
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Figure 9. Trajectories in the phase plane spanned by the maximal drop height and apparent 
contact angle for e = 10 -6 . Cases shown correspond to the ones in Fig.[8ja). 



A more complete picture of the time evolution for different initial profiles 
is obtained by considering the dependence of overall measures on time, and 
the trajectories of time evolutions in various "phase planes". For the axes 
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Time 

Figure 10. Shown is the dependence of the apparent contact angle on time for e = 10 -6 . The 
given cases correspond to the three evaporating drops in Fig.[8ja). 




1 2 3 4 1 2 3 



X X 

Figure 11. Given are for e — 10~ 6 drop profiles for selected drop volumes during the course 
of evaporation (for the three different initial profiles). We show as well the steady state drop 
of the same volume. Panel (a) gives the initial profiles at V = 1000 whereas panels (b) to (d) 
show profiles at V = 500, V = 100, and V = 10, respectively. 

of the latter we choose measures that do not change when the domain size 
is varied for an identical drop. We give results in two such phase planes, 
namely, the one spanned by volume and maximal drop height (Fig. [8]> and 
the one spanned by maximal drop height and apparent contact angle (Fig. [9]). 
The change of the contact angle over time is given in Fig.[l0| whereas Fig. \TT\ 



TTPll.tex; 11/04/2012; 2:46; p. 13 



14 



shows selected drop profiles. Figs.[8ja) and (b) compare results for very small 
e = 10 -6 and the largest used e = 1. As the results are qualitatively rather 
similar, the remaining figures [9]to [II] are for e = 1CT 6 only. 

Scrutinising Figs. [8] to [TT] one makes several observations: (i) The time 
evolutions starting from the three different initial profiles converge after some 
initial adjustments whose details depend on the particular initial profile shape, 
(ii) The convergence is slightly faster for smaller e. Here, "faster" means that 
the trajectories approach each other at higher volume [cf. Fig.(8j. In absolute 
terms the overall evolution becomes faster with increasing e. Thereby, the 
trajectories of the initial parabola and the initial profile taken from the steady 
state calculations approach each other earlier than they are approached by 
the trajectory of the initial gaussian. (iii) The family of steady profiles with 
influx represents drops clearly distinct from the "freely evaporating" shrink- 
ing drops. The family of steady profiles does not approach the trajectories 
of evaporating profiles when the drops become small. Even for very small 
drops their contact angle remains always larger by a roughly constant factor 
than the one of the evaporating drops. The factor is about two for e = 10~ 6 
and approaches four for e = 1. (iv) The overall picture in Fig.[8]for different e 
looks very similar, only the h max axes scale differently. A similar observation 
holds for the representations as given in Figs. [9] to 1 1 where, however, both 
axes would need to be scaled. 

Next, we discuss the behaviour of the different initial drop profiles at early 
times. We use Fig. [9] as an example. For larger e the behaviour is slightly less 
pronounced, but all the curves look qualitatively the same (not shown). For 
instance, one obtains a plot that is roughly the one for e = 1 when scaling the 
/i max -axis and # app -axis of Fig.[9]by factors of 1/6 and 1/40, respectively. A 
similar rule applies to Fig. [TUJ when additionally scaling time by about 1/5. 
In Figs.|9]and[T0| one finds for the initial parabola profile a strong decrease in 
the apparent contact angle at early times. This corresponds to an adjustment 
of the contact line region to the influences of the disjoining pressure. As the 
central drop region nearly coincides with the initial steady profile (per defi- 
nition at same volume and height) the two curves approach each other rather 
fast. In the course of the time evolution the central part of the profile remains 
a parabola. However, for the gaussian at early times the contact angle changes 
non-monotonically: The profile adjusts on the one hand its contact line region 
to the disjoining pressure influences (related to the "earlier wiggle" in the 
curve for the parabola in Fig. [9]). On the other hand, its central region adapts 
to a parabola (second "wiggle" in the curve in Fig. [9]). 

All three profiles approach each other after the initial adjustments. Their 
central part can be well fitted by a parabola, e.g., for V = 500 [V = 100] 
and e = 10~ 6 down to thicknesses of about h = 60 [h = 35]. Keeping 
the drop volume constant, that thickness decreases with increasing e and vice 
versa. In contrast, the steady profile of the drop with influx can be fitted by a 
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parabola in a smaller central part of the drop. The deviation from the parabola 
becomes clearly visible, e.g., for V = 500 [V = 100] and e = 10~ 6 at about 
h = 120 [h = 60] (already 20-30% below the maximum). This percentage 
range remains roughly the same when changing e for fixed drop volume. 

Our direct comparison of evaporating steady state drops with influx and 
evaporating drops without influx shows that the former cannot be used to 
approximate the latter as their shapes always differ. A freely evaporating 
shrinking drop has always a smaller apparent contact angle than the steady 
fed drop. This has been shown for a wide range of length scale ratios from 
e = 10~ 6 to e = 1. Note, however, that the differences slowly decrease for 
decreasing e. 




Figure 12. Space-time plots of an evaporating droplets for (a) e = 10 -6 , (b) e = 10~ 4 , (c) 
e = 10~ 2 , and e = 1.0. The initial profile is always the same parabola of maximal height 
-Hmax = 2500 and volume 4 x 10 5 on a precursor film of h = 1. The corresponding contact 
angle is ini = 20.8. 



We have observed that freely evaporating drops with a similar initial ge- 
ometry (volume and height) as per steady-state drops, for the range of e ex- 
plored, never spread macroscopically before their contact line recedes. How- 
ever, an initial spreading phase is often observed in experiments 151 l33l |4|. 
To investigate this further we perform a number of simulations starting with 



large parabolic drops. Fig. 12 gives a set of space time plots obtained for 



different length scale ratios from e = 10 to e = 1. All of them start 
from the identical initial profile. For small e < 10~ 3 [panels (a) and (b)] 
the behaviour is very similar to the one described above for drops with the 
same initial geometry as the steady drops: the drops shrink monotonically, 
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their height and width decrease slowly. However, at larger e > 10~ 3 [panels 
(c) and (d)] the behaviour is qualitatively different: At early times the drops 
spread. Thereby they loose height and gain width quite fast, the apparent 
contact angle decreases strongly. Then the drop reaches a maximal width 
before the contact line starts to recede again. In the shrinking stage, the height 
and width of the drop decrease slowly as before. The spreading is faster for 
larger e [Fig.fT^d)]. 

The contact angle for the initial profile is in all cases #i n i = 20.8. Com- 
paring this with Fig. |6jb) one notices that this angle roughly coincides with 
the limiting contact angle (for large drops) for e « 3 x 10~ 4 . This value lies 
between the regions (in e) where we find receding and spreading evaporating 
drops, respectively. Extrapolating from this finding, we formulate the hypoth- 
esis that the steady state drops with influx studied above in section [4] repre- 
sent limiting solutions between the case of spreading and shrinking freely 
evaporating drops (without influx). 




o 



8 



Figure 13. (colour online) Phase diagram indicating the initial behaviour of an evaporating 
drop. In the plane spanned by the initial contact angle 9^ and the length scale ratio e we 
indicate where the drop initially spreads, and where the contact line recedes right from the 
beginning. Each symbol corresponds to a time simulation. The solid line corresponds to the 
numeric result characterising the large steady drops with influx (cf. Fig.[6ftY)). 

To test the hypothesis we perform a number of time simulations with 
parabolic initial drops of different initial contact angle and at different e. All 
of them are of the same (large) volume. The results are given as a scatter 
plot in Fig. [T3] together with the curve for the large steady drops obtained in 
section [4] [solid line of Fig.^b)]. For each initial condition we record whether 
the drop spreads initially or directly starts to recede. Our results indicate that 
the above hypothesis seems to hold. The transition between initial spreading 
for large initial contact angle and a receding of the contact region right from 
the beginning roughly coincides with the power law dependence 6 ~ e -1 / 4 
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(Fig. |6jb), curve for large volume) not only in the power but as well in the 
prefactor. The prefactor of the curve obtained from the steady drops with 
influx seems to be slightly below the transition found in the time simulations. 
Further studies will be necessary to give a more detailed account. 



6. Conclusions 

We have analysed a thin film evolution equation for a wetting evaporating 
liquid on a smooth solid substrate. We have focused on slowly evaporating 
small sessile droplets where thermal effects are insignificant. Employing the 
model, we have first studied evaporating drops as steady state solutions for the 
case when they are fed through a porous part of the substrate. In particular, an 
asymptotic analysis has focused on the transition region between the precur- 
sor film and the bulk drop; and a numerical continuation of steady state drops 
has determined the fully non-linear drop profiles as a function of the overall 
influx for various values of the length scale ratio e. We have found that for 
large steady drops, the volume as well as the apparent contact angle decrease 
for increasing e roughly as e^ 1 / 4 . This does well agree with the scaling e -1 / 4 
determined via the asymptotic analysis. Note that the mentioned logarithmic 
corrections to both the overlap range where the matching is done and the 
resulting profile slope (apparent contact angle) are found to have opposite 
effects and are too subtle for order-of-magnitude estimates. 

Furthermore, we have employed the model to study the time evolution of 
freely evaporating drops that are not fed through the substrate, i.e., the full 
evolution equation has been numerically integrated. Thereby the time evolu- 
tion of several different initial drop shapes (for identical maximal height and 
volume) has been compared. It has been shown that freely evaporating drops 
with different initial profiles converge onto certain trajectories in phase space. 
However, a direct comparison of the freely evaporating drops without influx 
with the evaporating steady state drops with influx has shown that the latter 
cannot be used to approximate the former. A freely evaporating shrinking 
drop has always a smaller apparent contact angle than the steady drop with 
influx. Here it has been investigated for a wide range of length scale ratios 
from e = 10~ 6 toe = l. However, as the differences between the two types 
of profiles slowly decrease with decreasing e, further studies should scrutinise 
the case of even smaller e. 

We have noted that in our simulations the freely evaporating drops with 
a similar initial geometry (volume and height) to steady-state drops with in- 
flux, never spread macroscopically before their contact line starts to recede 
and the drop shrinks. As drops undergo an initial spreading phase in many 
experiments, we have investigated this further and found that drops spread 
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[shrink] from the beginning if their initial contact angle is larger [smaller] 
than the apparent contact angle of large evaporating drops with influx. 

This seems to be a very promising result that should be further scrutinised 
as it might have interesting consequences: (i) If the apparent contact angle 
of a steady drop with influx takes the role of an equilibrium contact angle 
9 e , relations between the dynamic angles and the contact line velocity known 
from non-volatile partially wetting liquids [6] could hold. Although, this has 
recently been shown for the case of evaporating partially wetting liquid [2], it 
remains an open question for the case of a wetting liquid that we study here, 
(ii) It might further be possible to predict the maximal drop radius and the 
contact angle at which the initial spreading ceases, and the 'turn around' to 
the receding motion occurs. It seems plausible that the profile at turnaround 
might actually be identical to the steady drop profile with influx at the same 
volume. Note that the freely evaporating drop spreads and evaporates, i.e., 
the volume at turn-around does not correspond to the initial one. 

Note finally, that here we have studied a particular evaporation model valid 
for small drops in situations where thermal aspects and the dynamics in the 
surrounding gas phase can be neglected. However, the non-isothermal models 
can all be studied with a similar methodology, i.e., the properties of steady 
drops with local influx can be determined and may be employed to gain a 
deeper understanding of the coupled transport and phase change processes. 
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